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REGULARITY OF SOLUTIONS TO STOCHASTIC VOLTERRA 
EQUATIONS WITH INFINITE DELAY 

ANNA KARCZEWSKA AND CARLOS LIZAMA 


Abstract. In this article we give necessary and sufficient conditions providing regularity 
of solutions to stochastic Volterra equations with infinite delay on a d-dimensional torus. 
The harmonic analysis techniques and stochastic integration in function spaces are used. 
The work applies to both the stochastic heat and wave equations. 


1. Introduction 

We consider the following integro-differential stochastic equation with infinite delay 

(1.1) X(t,6)= [ b{t-s)[AX(s,9) + ^^(s,9)}ds, t> 0, 9 G T d , 

J — oo OS 

where b G L * 1 (M_ ) _), A is the Laplace operator and T d is the d-dimensional torus. In 
OH), Wr is a spatially homogeneous Wiener process with the space covariance T and 
taking values in the space of tempered distributions S'(T d ). Such equation arises, in the 
deterministic case, in the study of heat flow in materials of fading memory type (see {^3, 

C33)- 

In this paper we address the following question: under what conditions on the covariance 
T the process X takes values in a Sobolev space H a (T d ), particularly in Ii 2 (T d )l 

This problem has been investigated, for the stochastic heat and wave equations in case of 
M 2 , by Dalang and Frangos j3] (see also Mueller USD- Their results were recently extended 
on W l and generalized to T d in [2]. In m, regularity of stochastic convolutions solving 
Volterra equations on M. d has been studied. The results have been obtained by using the 
resolvent operators corresponding to Volterra equations and expressed in terms of the 
spectral measure and the covariance kernel T of the Wiener process Wr- Our approach 
is different. We study a particular case of weak solutions under the basis of an explicit 
representation of the solution to (fOl (cf. formula EH). We remark that the knowledge 
of the regularity of solutions is important in the study of nonlinear stochastic equations 
(see e.g. @ and d). 

Observe that equation EU can be viewed as the limiting equation for the stochastic 
Volterra equation 

(1.2) X(t,9)= f b(t-s)[AX(s,9) + ^^(s,9)}ds, t> 0, 0 G T d , 

Jo OS 
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see mi Chapter III, Section 11.5] to obtain details on this assertion. If b is sufficiently 
regular, we get, by differentiating EH with respect to t, 


(1.3) ^(t, 9) = 6(0)[AX(f, 6) + ^(t, 9)} + jT b'(t - s)[AX(s, 9 ) 


dW r , mi , 

+ -^-(s,0)]ds, 


where t > 0 and 9 G T . 


Taking in (11.31) . 6(t) e 1 we obtain 
dX 

(1.4) l ~dt [t ' U) = ^^V) + '-dT 

\ X(0,9) = 0, 9eT d . 

Similarly, taking b(t) = t and differentiating (11.21) twice with respect to t we obtain 




(1.5) 


7 #’"> 

*( 0 , 0 ) 


= A X(t,0) + m^t,9), 


t > 0 ,9 e T d 


= 0 , 
= 0, 


9 e T d , 
0 G T d . 


It has been shown in |3, Theorem 5.1] (see also [HJ Theorem 1]) that equations (II .41) and 
(ESI) on the d-dimensional torus T d have an dd Q+1 (T (i )-valued solutions if and only if the 
Fourier coefficients (y n ) of the space covariance T of the process Wr satisfy 

(1.6) ^ 7 n (l + |n| 2 )" < oo. 

n£Z d 

Observe that for both, stochastic heat El and wave ESI equations, the conditions 
are exactly the same, despite of the different nature of the equations. On the other hand, 
the obtained characterization gives a natural framework in which nonlinear heat and wave 
equations can be studied. 

In this article, we will prove that condition EEI even characterizes H a+1 (T d ) - valued 
solutions for the stochastic Volterra equation EH), provided certain conditions on the 
kernel b are satisfied. This is a strong contrast with the deterministic case, where regularity 
of (11.11) is dependent on the kernel b. The conditions that we impose on b are satisfied by 
a large class of functions. Moreover, the important examples b(t ) = e -< and b{t) = 
are shown to satisfy our assumptions. 

We use, instead of resolvent families, a direct approach to the equation EH) finding an 
explicit expression for the solution in terms of the kernel b. This approach reduces the 
considered problem to questions in harmonic analysis and lead us with a complete answer. 


2. Preliminaries 

Let (Q, X, (Tt)t> o,-P) be a complete filtered probability space. By T d we denote the 
d-dimensional torus which can be identified with the product (— n, ir] d . Let D(T d ) and 
D’(T d ) denote, respectively, the topological vector space of functions (j) £ Cfi°(T d ) endowed 
with the usual topology corresponding to the convergence (see, {I or CH) and the space 
of distributions. By (£, </>) we denote the value of a distribution £ on a test function. 
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We assume that Wr is a D , (T rf )-valued spatially homogeneous Wiener process (abbr. 
s.h. Wiener process) with covariance F. For more information about s.h. Wiener process 
we refer to nni- 

Any arbitrary s.h. Wiener process Wr is uniquely determined by its covariance T E 
D'(t d ), which is a positive definite distribution, according to the formula 


(2.1) E(FF r (t, 9), 4>) (W v (t , 9), ip) = min(t, r) (r, cp * ip ( s) ), 

where <p,ip € D(T d ) and ip( s )(rj) = ip{—rj), for 77 E T d . Because Wr is spatially homoge¬ 
neous process, the distribution T = T(9 — 7 ) for 9,rj E T d . 

The space covariance T , like distribution in D'(T d ), may be uniquely expanded (see, 
e.g. 0 or [El ) into its Fourier series (with parameter w = 1 because the period is 2ir) 

(2.2) r(0) = ei{n,9 hn, 9 E T d , 


convergent in D'(T d ). In (12.21) . {n,9) = Yl < i=i n i9i and Z d denotes the product of integers. 

The coefficients y n , in the Fourier series El, satisfy (see, eg. Lesson 29, Section 
29.4] or 0 ): 

(1) = 7 _ n for n E Z d , 7 n are non-negative, 

( 2 ) the sequence ( 7 ^) is slowly increasing , that is 


(2.3) 



In 

1 + |n| r 


< + 00 , for some r > 0 . 


Let us introduce, by induction, the following set of indexes. Denote := N, the set 
of natural numbers and define Zf +1 := (Zj x Z d ) U {(0,n) : n E Z^}. Let us notice that 
Z d = Zf U (—Zf) U {0}. For instance, for d = 2, Z 2 S = N x Z U {(0, n) : n E Z}. 

Now, the s.h. Wiener process Wr corresponding to the covariance F given by El may 
be represented in the form 

(2.4) Writ, 9) = Snf3o(t) + Y [cos(n,9) (3i{t) + sin(n, 9) 0iit)\ , 


t > 0 and 9 E T d . 

In (12.41) Pq, pi, (3%, n E Z d , are independent real Brownian motions and 70 , 7 n are coeffi¬ 
cients of the series (12.21) . The series El is convergent in the sense of D'(T d ). 

Because any periodic distribution with positive period is a tempered distribution (see, 
e.g. 0), we may restrict our considerations to the space S'(T d ) of tempered distributions. 
By S(T d ) we denote the space of infinitely differentiable rapidly decreasing functions on 
the torus T d . 

In the paper we use the following definition of the Sobolev spaces. By H° = H a (T d ), 
a E R, we denote the real Sobolev space of order a on the torus T d . The norms in such 
spaces may be expressed in terms of the Fourier coefficients (see jl]) 

(2.5) lien*. = ( £ (l + |n| 2 )“|{„| 2 ) = (1 |{„| 2 + 2 £ (l + M 2 )» ((& + (£) 2 ) 

\n£Z d ) \ neZ d 



where = th + €=€- n , n E Z d . 
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There is another possibility to define the Sobolev spaces (see, e.g. d). 

3. Main results 

If 6 E T( oc (M + ) and fi E C, we shall denote by r(t,n) the unique solution in Lj oc ( R + ) 
to the linear Volterra equation 

(3.1) r(t, fi) = b(t) + n f b(t — s)r(s, n)ds, t> 0, 

Jo 

see [7J Theorem 3.1]. By taking Laplace transform in t we obtain r(X,fi) = 6(A)(1 — 
/x6(A)) -1 and hence, in many cases, the function r(f,/x) may be found explicitly. For 
example, if b(t) = 1, then r(A,/r) = or, inverting the Laplace transform, we obtain 
r(t,n) = e M *. Analogously, if b(t ) = e _t , then r(A,/i) = and hence r(t, n) = 

e (-i+/d*. Other examples are: 

sinh sjji t 


b(t) = t, 


b(t) = te _t , 

r (t, fi) 


sinh t 


d 7^ 0 


M 7^ 0. 


For more examples, see monograph m by Priiss. 

Let us denote by /(fc), k S Z, the fc-th Fourier coefficient of an integrable function /: 


f27T 


- 2n 


—ikt 




Given / defined on [0, 27 t] we denote again by / the periodic extension to R. Let b 6 
L 1 (M + ). We first observe from 

/ t POO 

b(t — s)f(s)ds = / b(s)f(t — s)ds 
-oo J 0 

that F is periodic of period T = 2ir as /. Now using Fubini’s theorem we obtain, 

(3.2) F(k) = b(ik)f(k), k € Z, 


where 


f c 

6 (A) := 

Jo 


e xt b(t)dt denotes the Laplace transform of 6 . 


In what follows we will assume that b(ik ) exists for all k E Z and suppose that A —> 6 (A) 
admits an analytical extension to a sector containing the imaginary axis, and still denote 
this extension by 6 . We introduce the following definition. 

Definition 3.1. Let r be given as in m- We say that a kernel 6 e L X (R + ) is admissible 
for equation eid ^ 


POO 

lim |n | 2 / [r(s,—\n\ 2 )] 2 ds =: Cb 

i|-*oo Jo 


\n 

exists. 

Example 3.2. i) In case b(t) = e~ l we obtain 

|2 


n 


POO 

/ [r(s,-\n\ 2 )] 2 ds = - 

Jo 1 


n 


+ \ n \ 
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and hence Cj, = 1. 


ii) In case b(t) = te t we obtain, after a calculation using e.g. Formula 2.662 (2)], 


n 


/»oo 

; / [r(s, — \n\ 2 )] 2 ds = 

Jo 


\n\ 


4 + 4|n | 2 


and hence Cb = 4 . 


By a solution X(t, 0) to the stochastic Volterra equation ( 11 . 11 ) we will understand a pro¬ 
cess X taking values in the space S'(T d ) and satisfying the equation (II . 1 II . The following 
is our main result. 


Theorem 3.3. Assume b £ L 1 (M + ) is admissible for d- Then, the equation d has 
an H a+1 (T d )-valued solution if and only if the Fourier coefficients (y n ) of the covariance 
T satisfy 

(3.3) ^ 7 n (l + |n| 2 )“ < oo. 

n£Z d 


In particular, equation d has an L 2 (T d )-valued solution if and only if 


E ln 

1 + Ini 2 

n&. d 


< oo. 


Proof. We shall use the representation (EH for the Wiener process Wr(t, 6) with respect 
to the basis: 1, cos(n, 0), sin(n, 9), where n £ T, d and 9 £ T d . Denote by W n (t,0) := 
cos(n, 9) Pn(t) + sin (n,6) /3 2 (t), n £ T , d , that is the n-th element in the expansion (12.411 . 
Equation Q may be solved coordinatewise as follows. 

Assume that 

(3.4) X(t, 9) = '^2[cos(n 1 9)X^(t)+sm(n,9)Xl(t)}-\-X 0 (t). 

n£ Z d 


Introducing El into El, we obtain 

cos(n, 9)Xn(t)+sin(n, 9)X 2 (t) = —\n\ 2 f b(t—s)[cos(n, 9)X^(s)+sm(n, 9)X 2 (s)]ds 

J — OO 

/ ‘lln [ b(t—s)[cos(n, 9) Pn( s )+ sin(n, 9) p 2 (s)\ds, 


or, equivalently 

[cos(n, 9), sin(n, 0 )] 


xk(t) 

X 2 (t) 




/ b(t — s)[cos(n, 9), sin(n, 0)] 

-OO 

t 


+ \/2~/ n / b(t - s)[cos(n,0),sin(n, 9)] 


^(»)1 ds 
* 2 (s)J ds 

’dPn(s) 

>2 


dPh( s ) 


Denoting 


x n (t) ■.= 


xk(t) 

Xfa) 


and Pn(t) : = 


Pk(t) 
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we arrive at the equation 

(3.5) X n (t) = —|n| 2 f b(t - s)X n (s)ds + yferfa f b(t - s)d(3 n {s). 

J — oo J—o o 

Taking Laplace transform in t. and making use of eu> with /j, = — |n| 2 and EHl . we get 
the following solution to the equation (13.511 : 


/ £ _ /*00 _ 

r{t - s.-ln^y/^r^dfdn^s) = / r{s,-\n\ 2 )^J 2 2q^d!3 n {t - s). 

-oo J 0 

Hence, we deduce the following explicit formula for the solution to the equation GU): 

_ no o 

X(t,e) = yfySPo(t) + cos (n,0) r(s,-\n\ 2 )d(3l{t - s) 

. J o 


(3.6) 




/o 

|2\ jo2, 


+ sin(n,0) / r(s, — |n| 2 )d(3 2 (t — s) . 

Jo 

From definition of the space H a and condition (12.511 the process X(t ) E iL a+1 , P-almost 
surely, if and only if 
(3.7) 

'/rt \' z / rt 

Ea + 


|R| 2 )" +1 7n 


(J r(t — s, —\n\ 2 )df3^(s)j + (J r(t - s,-|n| 2 )d/3 2 (s) 


<oo. 


In the representation n of Hr, processes PoiPniPni n ^ are independent real 
Brownian motions. So, the stochastic integrals in ETTll are independent Gaussian random 
variables, too. Hence, El holds P-almost surely if and only if 
(3 ' 8) 

J] (! + |n| 2 )“ +1 7nE[( [ r(t - s, -|n| 2 )d/^(s)) 2 + ( f r(t-s, -|n| 2 )d/3 2 (s)) 2 ] < c». 
neZ d J-oo J -oo 

Or equivalently, using L 2 isometry of stochastic integrals, if and only if 


(3.9) 


Ep + 


n| 2 r + 1 7n 


n oo 

/ [r(«,-H 2 )] 

Jo 


2 ds < oo. 


Since b is admissible for the equation o, we conclude that El holds if and only if 

^q+l 7rt 
Ini 


E t 1 + wT +I A < °°. 


and the proof is complete. 


Concerning uniqueness, we have the following result. 

Proposition 3.4. Assume b E L 1 (M + ) is admissible for El and the following conditions 
hold: 

(i) ^ 7n(l + |n| 2 )“ < oo, 


(ii) {l/b(ik)}kez CC \ {— |n| 2 : n E Z d }. 
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Then, El has a unique H a+1 (T d )-valued solution. 

Proof. 

Let X(t,9 ) be solution of 

X(t,6) = f b(t — s)AX(s, 0)ds. 

J — OO 

Taking Fourier transform in 0 and denoting by X n (t) the n-th Fourier coefficient of 
X(t,9)(t fixed), we obtain 

X n (t) = — |n| 2 f b(t — s)X n (s)ds 

J — OO 

for all ngZf Taking now Laplace transform in t, we get that the Fourier coefficients of 
X n {t ) ( n fixed) satisfy 

(l + |n| 2 &(ifc))X n (A;) = 0 

for all fe G Z. According to (ii) we obtain that X n (k) = 0 for all k G Z and all n G Z rf . 
Hence, the assertion follows by uniqueness of the Fourier transform. 


The following corollaries are an immediate consequence of Theorem Id. Ml The arguments 
are the same as in |5i. We give here the proof for the sake of completeness. 

Corollary 3.5. Suppose b G L 1 (R + ) is admissible for El and assume T G L 2 (T d ). Then 
the integro-differential stochastic equation El has a solution with values in L 2 (T d ) for 
d = 1,2,3. 

Proof. We have to check equation (13.31) with a = — 1. Note, that if T G L 2 (T d ) then 
f = (j n ) G l 2 (Z d ). Consequently 

^ 1 + |n| 2 - ^ 7n 2 ) 1/2 (J2 (l + |n| 2 ) 2 ) 17 "' 

n£Z d n£Z d 

But I2nez. d Tn 2 < oo and 

(3.10) 7 -i— tkt- < oo if and only if 2 q > d. 

K ^ 1 + \n 2 )i J H 

n€ z d 

Hence, the result follows. 


Corollary 3.6. Suppose b G L 1 (R + ) is admissible for El and assume f G l p (7j d ) for 
1 < p < 2. T/ien t/ie integro-differential stochastic equation El has a solution with values 


in L 2 (T d ) for all d < ^EL. 
Proof. Note that 


E 

n£ Z d 


7n 


1 + In 




(E E 


(1 + |n| 2 ) < ? 


) 1/9 , 


where - + - = 1. Hence the result follows from (ETUI) with q = -£[■ 


For a = — 1, the condition USD can be written as follows. 
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Theorem 3.7. Let b be admissible for the equation m. Assume that the covariance T 
is not only a positive definite distribution but is also a non-negative measure. Then the 
equation tn has L 2 (T d )-valued solution if and only if 


(3.11) 


(r, Gd ) < +oo, 


where 


(3.12) 



x E T d . 


The proof of Theorem 13.71 is the same that for Theorem 2, part 2) in jH|, so we omit it. 
The functions Gd(x ) are apparently well-known. If d = 1, then G\ is a positive con¬ 
tinuous function. If d = 2, G% is continuous outside 0 with singularity at 0 of the form 

c log |—r, c > 0. For d > 2, Gd(x) ~ c,—-pvr as x — * 0. For more details concerning the 

3/ v / \X ^ 

function Gd we refer to 0 Prop. 3] and [TT|- 

From the above properties of functions Gd defined by (ET21) and the condition G3D 
we obtain the following result (see |5, Theorem 4]). 


Corollary 3.8. Assume that T is a non-negative measure and b is admissible. Then 
equation EJ has function valued solutions: 

i) for all T if d = 1; 

ii) for exactly those T for which J^| <:L In |0| T(d0) < +oo if d = 2/ 

in) for exactly those T for which jj 6 ,| <1 T(d0) < +oo if d > 3. 


In what follows, we will see that formula m also implies that X is Holder continu¬ 
ous with respect to t. In order to do that, we need assumptions very similar to those in {2^. 

Hypothesis (H) 

Assume that there exist 5 E (0,1) and Cs > 0 such that, for all s E (—oo, t) we have: 


(i) 


[r(t - t,-\ n\ 2 )] 2 dr < C 5 \n\ 2 ^ \t - s\ 5 ; 


•s a S 

(ii) f [r{t — t, — \n\ 2 ) — r(s — r, — |n| 2 )] 2 dr < \t — s\ 

J — OO 


Proposition 3.9. Assume that — 


7 n 




+ \ n \ 


< +oo. Under Hypothesis (H), the trajec¬ 


tories of the solution X to the equation E2D are almost surely p-Holder continuous with 
respect to t, for every g E (0, 5/2). 
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Proof. From the expansion EH1 and properties of stochastic integral, we have 
E \\X(t, 0 ) - X(s,0)\\l 2 = E || v /7o(/3 0 (t) - /3o(s)) 

+ c °s(n, 6) f f r(t-T,-\n\ 2 )d/3l(T)- f r(s-T,-|n| 2 )d^(T) 

L \J ~ OO J— OO 


L 2 


+ sin (n,0)l r(t—T,—\n\ )df3 n (r) — / r(s-r,-|n| )d/? ri (r) 

= (2vr) d (701i - s| 

+ ^ 27 n f [r(t— r,— |n| 2 ) — r(s — r,—| \n\ 2 )] 2 dT + f r 2 (t—T,—\n\ 2 )dT 
n£Z d U-00 Js 

According to assumptions (i) and (ii) of the Hypothesis (H), we have 

E\\X(t,0) - X(s, 0)\\ 2 l2 <C s J2 2 ln \n\ 2 ^\t-s\ 5 . 

nGZf 

Because A is a Gaussian process, then for any m £ N, there exists a constant C m > 0 
that 


E\\X(t,d)-X(s,0)\\ 2 L 7 <C, 


Cs ^2 


n 


2(5-1) 




Taking m such that m5 > 1 and using the Kolmogorov’s criterion for continuity (see, e.g. 
EB3 ), we obtain that the solution X(t,6) is ?)-Holder continuous, with respect to t, for 
r] = 5/2 — 1/(2 m). ■ 


Example 3.10. Let us consider the particular case b(t) = e t > 0. Then, by previous 
considerations, r(t, — \n\ 2 ) = e( -1— l n l ' )t . One can check that in this case the Hypothesis 
(H) is fulfilled. 

Remark 3.11. We observe that the condition (i) in Hypothesis (H) is the same as 

\n\ 2 f [r(s, — |n| 2 )] 2 ds < C$\n\ 28 \t\ 5 

Jo 

and hence it is nearly equivalent to say that the function b is admissible. 

Acknowledgement The authors would like to thank the anonymous referee for the 
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